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Abstract 

We investigate a flexible method in which we can test the unitarity of 
the quark flavor mixing matrix step-by-step. Singular-Value-Decomposition 
(SVD) techniques are used in analyzing the mixing matrix over a broader 
parameter region than the unitarity region. Unitary constraints make us ex¬ 
tract CP violating properties without any specific parametrization when the 
magnitudes of at least three mixing matrix elements in three generation quark 
mixing are given. This method can also be applied to the analysis of lepton 
flavor mixing, in which only a few moduli are presently measured. 
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I. INTRODUCTION 


The Cabibbo-Kobayashi-Maskawa (CKM) matrix makes us possible to explain all 
flavor changing weak decay processes and CP violating phenomena up to now. Unitarity of 
the CKM matrix in the standard model (SM) is a unique property that we cannot loosen. 
We can use any parametrization of the CKM matrix as long as its unitarity is conserved. 
The original parametrization for three generation quark mixing is the Kobayash-Maskawa 
(KM) parametrization. The standard parametrization proposed by Chau and Keung 
is the product of three complex rotation matrices which are characterized by the three 
Euler angles, 6*12, 6^13, 6*23 and an CP-violating phase ^13. More widely used one is the 
Wolfenstein parametrization p, which was suggested as a simple expansion of the CKM 
matrix in terms of the four parameters: A, A, p and p. It has been also known that the 
CKM matrix for the three-generation case can be parameterized in terms of the moduli of 
four of its elements 0 . This four-value-KM ( 4 VKM) parametrization is rephasing invariant 
and directly related to the measured quantities. In three generation case we always need 
four independent parameters to define a unitary 3 x 3 matrix, as explained, eg. 6*12, 6*13, 6*23 
and 5 i 3, or A, A, p and p or even only moduli of any four independent elements of the 
matrix. 

The 4 VKM parametrization has several advantages over the other parametrization. This 
parametrization doesn’t need any specific representations for the mixing angles as long as 
the CKM is unitary, and no ambiguity over the definition of its complex CP phase is present 
above all. Secondly, the Jarlskog invariant quantity Jcp and non-trivial higher invariants 
can be reformulated as functions of moduli and quadri-products 0 . However, in the 4 VKM 
parametrization initial four-moduli input values should be fixed by experiments. Once we set 
four moduli to specific values, remaining five moduli of mixing elements are automatically 
fixed and we may lose some characteristic effects from interplaying between the moduli. In a 
conceptual point of view it is better if we can reduce number of a priori experimental input 


2 


values. This paper presents a novel parametrization in which we start with three-moduli 
input values. Through simple algebraic relations we can determine remaining six moduli 
of mixing elements. With more broader parameter space we can check the compatibility 
between measured values of mixing elements and their unitarity properties step-by-step. 

Many groups have made global hts and numerical works on CKM matrix elements with 
conventional representations which satisfy unitarity |^. One of the problems in these con¬ 
ventional parameterizations is that they are fully and completely unitary and are not flexible 
to include possible non-unitary properties resulted from unknown new physics. Therefore, it 
is a complicate task to make a step-by-step test to check the unitarity with experimental data 
if you use a unitary parametrization. In the following, we present three extended dehnitions 
for the unitarities of mixing matrix V in the order of the strength of the constraints: 

• Weak Unitary Conditions (WUC): We dehne that the mixing matrix V is weak uni¬ 
tary if it satishes 

\Vicf\^ = 1 i = u,c, t, and (3 = d,s, b. (1) 

" j 

These constraints appear to be well satished experimentally for the three generation 
case, and we start from this. Actually it was pointed out that there is an apparent 
functional violation in the available data: |I4dP + < 1 |0- For such a 

case with J 2 a = a < 1, we can easily generalize our method, and we just start 

with this new condition. 

• Almost Unitary Conditions (AUC): In addition to the constraint Eq. (P, if the fol¬ 
lowing constraints are satished 

Vi*Vja = = 0 for some parts of i, j = n, c, t, and a, f 3 = d, s, b, (2) 

let us call the mixing matrix almost unitary. Some combinations, which do not satisfy 
Eq. ( 2 ), may not make closed triangles, and may have different areas even though 
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making closed triangles. We have no specific models in which the mixing matrix 
satisfies this almost unitary conditions. Therefore, we will not consider the case with 
AUC. 

• Full Unitary Conditions (FUC): This corresponds to usual unitarity in which 
Eqs. ( 0 )> ( 0 ) are satisfied for all the indices. All six unitarity triangles from Eq. 
have the same areas. 

In Sec. II, we propose an alternative and more flexible parametrization of the CKM ma¬ 
trix in terms of the three moduli and the one independent parameter, which is induced by the 
singular-value-decomposition (SVD) method. We describe how to get the new parametriza¬ 
tion of the CKM matrix by using the SVD method in the three-generation case. Unlike 
from the previous parametrization with four moduli |^, we have more flexible leverage to 
test the unitarity step-by-step. We start with only three moduli rather than four moduli, 
and the remaining one can be adjusted depending on the condition of the unitarity, which 
we apply, i.e. WUC or FUC. In Sec. Ill, we analyze the CKM matrix numerically with our 
parametrization with the SVD method. Conclusions are also in Sec. III. Appendices A-B 
include details about the SVD method. 


II. NEW PARAMETRIZATION OF THE CKM MATRIX BY THE SVD METHOD 


We start with a definition in such a way that it satisfies the weak unitary conditions, 
Eq. (1): we have six constraint equations for the 3 generation mixing. These constraints 
are only parts of unitarity conditions and the introduced mixing matrix V may not be fully 
unitary. We study this explicitly with three generation quark flavor mixing matrix V in their 
absolute values and choose three independent moduli as starting input parameters. Explicit 
analysis depends on the choice of three input parameters. We consider the case with: 
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• Our Choice (Set-A): Input parameters |I4s|, \Vub\) \Vcb\- 

We can also choose different sets of input parameters, as examples: 

• Set-B: Input parameters |I4d|, |Ks|, |Wd|, 

• Set-C: Input parameters |I4s|, iWsI, \Vcb\- 

Mathematically three parametrizations of Set-A,B,C are all equivalent if three input values 
of each set are independent one another and all equally precisely measured. However, in 
reality, the upper-left 2x2 part of CKM matrix is approximately unitary and only one 
independent variable is dominantly evident, for example, the parameter A in the Wolfenstein 
parametrization or the Cabibbo angle sin 0c- Therefore, Set-B would be the worst choice 
for numerical analyses. For our choice of Set-A, the three inputs are all off-diagonal and 
independent each other, and all three values can be determined by three semileptonic decays, 
in which new physics contributions are severely suppressed. Therefore, we select upper off- 
diagonal elements in V, namely, |I4s|, iWfel, |Wfe| as the initial input variables in our analysis, 
i.e. the case with Set-A. 

If we are given the three input values of Set-A, then we get the values \Vud\ and \Vtb\'- 

\Vud? = l-\Vus?-\Vub\\ (3) 

\Vtb? = l-\Vub?-Kb?. (4) 

To obtain four remaining elements,|I4rfI, iWsI, \ Vtd\ and 114^1, we write four constraints for 
these four elements in Eq. (P as a matrix form: 

RX = B, (5) 


where 
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^ 1100 ^ 

0 0 11 

R = 

10 10 
^ 0101 ^ 

x = {\vj^,\Vcs\MVtd\MVts\Y, 

B = {1- iKbP, 1 - \Vtb\\ 1 - \Vud\\ 1 - \Vus\y. 


( 6 ) 

(7) 

( 8 ) 


In Eqs. (^, (P), X and B are column vectors and T means transpose of the matrix. Because 
of deti?=0, there is not a unique solution if any. In such a situation there exists a very 
powerful set of technique, known as Singular Value Decomposition (SVD) method. The 
details of the method are given in the Appendix A. Remaining mixing elements are expressed 
as follows: 

|ld:d| ® T Ui, 

|hcs|^ = a + U2, (9) 

= a + M 3 , 

+ It-A, 

where 

Mi = ^(i+ 2 |K.r+iK 6 r- 2 iKb| 2 ), 

U2 = ^(3 - 2|K/ - iK^r - 2|IA6|'), (10) 

U 3 = + 2|I4is|^ + 3|I4ife|^ + 2|V);fe|^), 

M4 = ^(i- 2 iK/+iKbr+ 2 iK 6 | 2 ), 

and a new variable, ‘a’, is introduced as a coefficient attached to the general solution. If 
there is no flavor mixing, we can set a = 1/4. The value of ‘a’ can be determined from 
Eq. (9) if we know any one value of \Vcd\i |h/s|, iT/dl, |h)s|- Constraints of non-negative |VjP 
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are applied for the range of variable a: 


^min niax( U2i ^max min^Mi,'U 4 ). 

We note that when three inpnt valnes |14s|, |14fe|, \ Vcb\ are given, the modnli sqnared of re¬ 
maining fonr mixing elements \ Vcd\, |Ws|, \ Vts\ are jnst qnadratic fnnctions of parameter 

a. As ‘a’ increases, |\4s| and \Vtd\ increase while |V"cd| and |14s| decrease. \Vud\ and \Vtb\ are 
hxed by the three inpnt values and are independent from the parameter a. And the bounds 
on the parameter ‘a’ will determine the regions of FUC and WUC, which will be explained 
later. 

As a next step, we further assume that the mixing matrix V satishes full unitary condi¬ 
tions. Then we have six more constraints: 

H = 0, (i, k) = (n, c), (u, t), (c, t), 

j=d,s,b 

E = {t,k) = {d,s),{d,b),{s,b). (11) 

j=u,c,t 

These constraints cannot be represented without introduction of complex numbers analyt¬ 
ically. If we know all the absolute values of V, however, we can express necessary and 
sufficient conditions for the constraints, Eqs. (0), in a geometric way. Eqs. (0) give six 
unitarity triangles corresponding to each six constraints, and all six triangles have equal 
area that is directly related to the Jarlskog’s rephasing invariant parameter Jcp- If we take 
one of the constraints Eqs. 0, for example, 

j=u,c,t 

a triangle is composed of three sides with lengths |14d||I4fe|, |Vcd||K:fe|, and |V)d||I4b|, with a 
necessary condition 


WdllWfel < \Vud\\Vub\ + \Vtd\\Vtb\, 


( 12 ) 
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where the equality holds in CP conserving case. For more general argument, let us rewrite 
Eq. (|^ as follows: 


^2 < ^1 + h , ( 13 ) 

where, as an example, h = | 14 d||KfeM 2 = I Kd 11 Kb I, and I 3 = |Kd||Kb|- After taking the 
square on both sides of Eq. (|I^ we can rearrange the constraint equation as follows: 

fih, h, h) = 2llll + 2llll + 2llll -It- It - It > 0, (14) 


where we denote newly introduced function / for later use. Using the Heron’s formula, the 
square of triangular area can be rewritten as follows: 


— s{s — h){s — l 2 ){s — h) — T7/(^i! ^ 2 , h)y 

lb 


(15) 


where s = {li + I 2 + h)/2. So the necessary condition ([T^) for the complete triangle means 
non-negative value of A^. The Jarlskog’s invariant parameter is written as follows: 


Jcp — 2 A— -\Jf{h, I 2 , h)- 


(16) 


If we expand / in terms of parameter a: 

/ = -(l-|KbP)V 

+2[Kd\^Kb\\\Vtt\^ - |Kb|2)(|Kfe| V - |K,|K3)(|Kfe|2 + \Vth\^)]a 

+2\Vud\^\Vub\^[\Vcb\\i + \Vtb\\3] - (iKfelV - \Vtb\\3f - |Kd|1Kb|^ (17) 


where the function / is quadratic of a. We can get the boundaries of the constraint Eq. o. 
and denote two roots of the quadratic equation as a_ and a+(> a_). Two real roots and 
the boundary points in the interval depend on only the three input values. Nonexistence of 
real solutions of the quadratic equation means that the three input values do not allow the 
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FUC. We note that if we force the mixing matrix V to be fully unitary, then six triangles 
from the constraints Eq. (ph]) have the same area, which are the sufficient conditions for the 
FUC. 


We can relate the coefficient a to the CP violating parameter in another representation 
of mixing matrix with the FUC. Let us consider the standard parametrization of the CKM 
matrix 


VcKM — 


C 12 C 13 


S 12 C 13 


sise 


-i(Si3 


-S12C23 - Ci2S23Sl3e*^^® C12C23 - Sl2S23Sl3e*'^^® S23C13 

S12S23 - Ci2C23Sl3e*'^i® -C12S23 - Si2C23Si3e*'^i3 C23C13 


(18) 


where Sij = sm.9ij,Cij = cos%. We hnd that the coefficient a is directly related to the 
parameters in the standard representation: 

a = -2Si2Ci2S23C23Si3COs5i3 - COS 2012 COS 2023 [“3 + COS 2013]. (19) 

o 


Three angles a, (3 ,7 of the unitarity triangle, which characterize CP violation, are dehned 
as follows: 


a = Arg[-(VuV,l)/{V^,V:t)], ( 20 ) 

13 = Arg[-(V,J^X)nVuVi)\, ( 21 ) 

7 = Ara[-(V^V:^)HVJA:Ci]. ( 22 ) 


The sum of those three angles, dehned as the intersections of three lines, would be always 
equal to 180*^, even though the three lines may not be closed to make a triangle, i.e. in case 
that CKM matrix is not unitary at all. We can also dehne these quantities from the area of 
the unitary triangle and its sides: 


sin/3' 


2 A 

\Vu\\Vtb\\V,d\\Vci> 


( 23 ) 
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(24) 

(26) 


sin 7 

a 


/ 


/ 


2A 

~ \Vud\\Vub\\Vj\Vcb 
= TT - - y, 


when the FUC is fully satished and the area of the triangles can be dehned from ([T5|). Any 
experimental data that indicates a y a' or y /?' or 7 y 7 ' means that three generation 
quark mixing matrix V is not fully unitary. 


III. NUMERICAL RESULTS AND DISCUSSIONS 


For given input values of |14s|, \ Vub\^ |fcfe|, the parameter ‘a’ is divided into two regions 
depending on whether the FUC is satished or not. We can divide the range of the parameter 
into two by setting h = \Vud\\Vub\,l 2 = iWdllWfeMs = as an example: 

• Region I : The maximum among li,l 2 ,h is larger than sum of the other two values. 
In other words it is not possible to make any triangle with these three segments. This 
region is outside of the interval of (a_, a+). 

• Region II : The maximum among li, ^ 2 , h is smaller than sum of the other two values. 
In other words it is possible to make a unitarity triangle. This region is conhned to 
(a_,a+). 


In region I, we cannot dehne Jcp- On the contrary we can dehne Jcp in region II and 


calculate it with Zi, Z 2 , h as shown in Eq. (^). In general the region II is surrounded by the 
region I. Two boundary points of region II correspond to the case of CP conserving case. 


For numerical analyses, we refer to the Particle Data Group (PDG) Gurrent values of 
three input moduli and corresponding sources of measured matrix elements are summarized 
in Table I. The input values of |Uus|, iWfel, \Vcb\ are randomly generated within 95% GL with 
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TABLE I. Input values of the matrix elements and their sources referred from the PDG. The 


output values are allowed intervals (95% CL) for WUC and FUC. 



matrix elements 

PDG values 

Sources 




iKsI 

0.2196 ±0.0026 

Kes decays 



Input 

\Vub\ 

(3.6 ±0.7) X 10-3 

B semileptonic decays 




\Vcb\ 

(41.2 ±2.0) X 10-3 

B semileptonic decays 




matrix elements 

WUC 

FUC 

PDG 



\Vcd\ 

0.210 ~ 0.224 

0.214 ~ 0.224 

0.219 ~ 

0.226 


IVcsl 

0.9735 ~ 0.9768 

0.9735 ~ 0.9760 

0.9732 ~ 0.9748 

Output 

\Vtd\ 

0.004 ~ 0.045 

0.004 ~ 0.014 

0.004 ~ 

0.014 


\yts\ 

0.001 ~ 0.045 

0.035 ~ 0.045 

0.037 ~ 

0.044 


uniform distributions. Each input determines both two regions for the WUC and FUC. 
The WUC is conhned to the interval (omm, cimax) of which calculations are described in 
the previous section. The FUC is conhned into the interval (a_,a+) which is a subset of 
{cimin, dmax)- In the restricted regions we again generated randomly the valnes of parameter a 
for onr nnmeric calculations. Fig. 1(a) shows scattered points for \Vtd\ and |Lts| valnes when 
we apply the WUC to the choice of parameter a. The scattered points compose a quadrant 
in the lUdhIhtsI plane. In the hgure we also draw the curved axis for the parameter a. The 
labels on the cnrve are valid only when we set the three inputs to the center values in Table 
I. If other input values are taken, the numeric labelling should be slightly changed. Fig. 
1(b) presents scattered points when we apply the FUC to the choice of parameter a. The 
allowed region for FUC is mnch narrower than that for the WUC and is included in the 
region for the WUC. The curved axis for the parameter a is identical to that in Fig. 1(a). 
Fig. 2 shows corresponding scattered points for \Vud\ and \Vcd\ when we take the WUC and 
FUC. The moduli \Vud\ does not depend on the parameter a and the \Vcd\ is directly related 
to the valne of a. The axis for the parameter a is, therefore, a vertical line along the axis of 


11 












FIG. 1. The scattered plots for \Vtd\ and 114^1 values which satisfy the WUC (a) and FUC (b) 
in the case of Set-A. The allowed points are calculated from uniformly generated three input values 
in 95% CL. The curved axis for the parameter ‘a’ is drawn with the centered input values in Table 
I, not with the randomly generated input values. 
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IKdl. 

In Table I, we show the numerical output values for the moduli, \Vcd\) |Ks|, |hts|, 

within 95% CL with the unitary conditions, WUC and FUC. For comparison, we also show 
the results of the PDG values. The PDG values of individual matrix elements were deter¬ 
mined from three-level constraints from weak decays of the relevant quarks or from deep 
inelastic neutrino scattering together with the assumptions of 3 generation FUC. As can be 
seen, the allowed regions for the WUC are much broader than those for the FUC and the 
latter are subsets of the former. Our numerical results are consistent with 90% CL on the 
magnitudes of mixing elements in PDG: \Vtd\ — 0.004 ~ 0.014 and |Vts| — 0.037 ~ 0.044, 
particularly. We can see that the FUC and hierarchical input values of iWsI, iWbl, iWftI imply 
|htd| < lUsI- This contrasts with the results of the WUC. In this case, it is possible that 
\Vtd\ is equal to or even larger than \Vts\- 

If we start with different mixing elements, like |I4d|, |Ws|, |Wd| (he. Set-B), then \Vtd\ 
and \Vuh\ are hrst hxed and we introduce new parameter b from the SVD method, as shown 
in Appendix B in detail. Remaining four moduli, |Ucs|, |h7s|, and |Ub| are dependent on 
the parameter b and have correlated values one another. Precise measurements for one of the 
four moduli will £x remaining three moduli. Similarly, if we start with |Uus|, |Ucs|, iWbl (he. 
Set-C), then \ Vcd\ and \ Vts\ are hxed and new parameter c will be introduced by following the 
SVD method. Remaining four moduli, |I4d|, |Wb|, 1^7^! and |Ub| depend on the parameter 
c and have correlated values one another. However, as explained in Section II, it would 
be much more difficult numerically to analyze with Set-B or C compared to Set-A due to 
approximate unitarity of the upper-left 2x2 part of CKM matrix. 

To conclude, we proposed a hexible method in which the unitarity of quark-mixing 
matrix can be tested step-by-step. The singular-value-decomposition (SVD) method is used 
in analyzing the mixing matrix over a broader parameter space than the unitary region as 
well as in presenting a new parametrization of the CKM matrix. The question whether the 
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mixing matrix satisfies the WUC or FUC is a quite difficult and complicate matter within the 
standard PDG parametrization or similar unitary parametrization. In the parametrization 
by the SVD method the CKM matrix is represented with three moduli and an additionally 
induced flexible parameter 00. Once the value of the induced parameter a is determined, we 
can easily distinguish the FUC from the WUC. For example, with the Set-A input data, 
if we get 0.224 < a < 0.226, the mixing matrix satishes the FUC within 95% CL. From 
Fig. 1(b) we can also conclude that the FUC violated if \Vtd\ > 0.02 or \Vtd\ < 0.004. If 
\Vtd\ > 0.05, even the WUC is not satished. Fig. 2(b) shows that there is a strong correlation 
between \Vud\ and \Vcd\- This method can also be applied to the analysis of lepton flavor 
mixing, in which only a few moduli are presently measured. 
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Appendix A. THE SINGULAR-VALUE-DECOMPOSITION METHOD 


For detailed description of this method we consider the specihc case with the inpnt 
parameters, |14s|, |Uc6|, which is in Sec. II. In this case we have to solve Eq. (i- 

According to the method, the matrix R can be decomposed as a prodnct of three matrices: 

R = UWV^, (A.l) 


where 


U 


W 


V 


( 


1 

2 


v 


-11 1 -1^ 
-1 -1 -1 -1 
-11-11 
- 1-11 1 ^ 


^20 00 ^ 
0 72 0 0 
0 0 72 0 
^00 00 ^ 




/ 


1 

2 


V 


-1 72 
-1 0 
-1 0 
-1 -72 


72 -1 
-72 -1 


(A.2) 


(A.3) 


(A.4) 


In general the matrices IJ and V are orthogonal in the sense that their colnmns are orthonor¬ 
mal. 


^ V UikUiyi ^kW) 

(A.5) 

^ ^ ^ik^in ^kw 

(A.6) 


i 


We note that this decomposition is not nnique. For further practical calculations of SVD we 
refer to Ref. [^. The solutions of Eq. (|^) are obtained in two types, special solution and 
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general solution, which can get in two different ways. First, the special solution according 
to SVD is calculated by defining the inverse of R as follows: 


R = V[diag{l/wii)]U^ 

^3 -1 3 -1 ^ 
3-1-13 
-13 3-1 

- 13-13 


(A.7) 


(A.8) 


where we take l/wa = 0 if wu = 0. The matrix R is unique and does not depend on the 
way how the matrix R is decomposed. The special solution Xg is 


X, = RB. 


(A.9) 


The inverse matrix R does not satisfy the constraints which must be obeyed in the general 
sense of the inverse, namely, 

RR^RR^ I. 


However it satisfies 


RRB = B. 


Therefore, we can introduce general solutions Xg such that 

RXg = 0, (A. 10) 

and we can add this to the special solution Xg. We can see the general solutions of Eq. (|A-10|) 
by simple guessing as follows: 

X, = a(-l,l,l,-l)^ (A.ll) 
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where coefficient a can take any real value. Actually the coefficient a must be further 
conhned in such a way that the values of the mixing elements Vij should be within the range 
[0,1]. We can express the complete solutions as 

X = Xg + Xs. (A. 12) 


In algebraic terms, Eq. (^) dehnes i? as a linear mapping from the vector space of X to the 
vector space of B. If R is singular, then there is some subspace of X, called the nullspace, 
that is mapped to zeTo,RX = 0. The number of linearly independent vectors that can 
be found in Eq. (|A-10|) is the dimension of the nullspace called the nullity of R. In three 
generation quark mixing case the nullity is 1. 


Appendix B. CASE WITH SET-B INPUT PARAMETERS 


If we take three independent input parameters as \Vud\, |I4s|, and |Wd|, we can apply 


the same procedure to obtain remaining mixing elements. In this case the matrix Eq. 
becomes RX = B with 


^ 1010 ^ 

0 10 1 

R = 

110 0 
^0 0 11 ^ 

X = {\v,,\Mv,b\MVts\MVtb\Y, 

B = {l- IK.P, 1 - iKfeP, 1 - iWdP, 1 - \Vtd\Y. 


(B.13) 


(B.14) 

(B.15) 


Two obvious relations from the weak unitary conditions are: 

\vtd\^ = 1 - iK^r - iK^r. 
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(B.16) 

(B.17) 







Following the procedure described in the Appendix A, we can write total solution of X as 
follows: 


\Vcs\^ = b + wi, (B.18) 

\Vcb\'^ =—b-\-W 2 , (B.19) 

\Vts\‘^ =—bW 3 , (B.20) 

= bW 4 ^, (B.21) 

where 

Wi = ^(3 - iK^r - 2|K.P - 2|Kdr), (B.22) 

W 2 = ^(1+iK^r+ 2 iK.r - 2 iKdr), (B.23) 

W3 = ^(1 + iKdP - 2|K.P + 2|14rf|2), (B.24) 

w, = ^(-1 + 3|KdP + 2|K.P + 2\VJ^), (B.25) 


and the b is newly introduced parameter. In this case, if there is no flavor mixing, we can 
set b = 1/2. Eq. (|B.21|) shows that any values of |V/s|, |V/b|, \ Vts\-i \ Vtb\ will determine the 
value b. Further constraints are applied for the range of parameter 6 by > 0: 6mm = 
max{—Wi, — W4), bmax = min{w 2 , W3). For the FUC we expand / in terms of parameter b: 


/ = -(!- Kdl^b^ 

+2[\Vud\"Kb\\\Vt,\^ -\VJ^){\VJ^W2 - V)(|Kdr + \Vtd\")]b 

+2\Vudf\Vubf[\V,dfw2 + \Vufw4 - {\V,dfw 2 -\Vtdfw 4 f -\Vudf\Vubf. (B.26) 


Like the previous case we can get the boundaries of the constraint Eq. (p!3[) , and denote two 
roots of the quadratic equation as 6_ and 6+(> 6_). 
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